A nonextensive entropy approach to kappa-distributions 



Manfred P. Leubner 

Institute for Theoretical Physics, University of Innsbruck 
A-6020 Innsbruck, Austria 



Most astrophysical plasmas are observed to have velocity distribution func- 
tions exhibiting non-Maxwellian suprathermal tails. The high energy particle 
populations are accurately represented by the family of kappa-distributions 
where the use of these distributions has been unjustly criticized because of 
a perceived lack of theoretical justification. We show that distributions very 
close to kappa-distributions are a consequence of the generalized entropy fa- 
vored by nonextensive statistics, which provides the missing link for power-law 
models of non-thermal features from fundamental physics. With regard to the 
the physical basis supplied by the Tsallis nonextensive entropy formalism we 
propose that this slightly modified functional form, qualitatively similar to the 
traditional kappa-distribution, be used in fitting particle spectra in the future. 



1 Introduction 

A variety of space observations indicate clearly the ubiquitous presence of suprathermal 
particle populations in astropliysical plasma environments 0. The family of kappa ve- 
locity space distributions, introduced first by Vasyliunas 0, is recognized to be highly 
appropriate for modeling specific electron and ion components of different plasma states. 

Numerous magnetospheric interaction processes and instabilities were studied success- 
fully within the concept of kappa-distributions 0, Q ranging from plasma sheet ion and 
electron spectra ^ to a combined electron-proton-hydrogen atom aurora The sat- 
uration of ring current particles towards a kappa-distribution was studied by Lui and 
Rostoker and accurate fittings of observed electron flux spectra were performed by a 
power law at high energies [Q. Furthermore, it was demonstrated by Leubner that 
the Jovian banded whistler mode emission can be interpreted within a kappa-distribution 
approach and that mirror instability thresholds are drastically reduced in suprathermal 



space plasmas [0, |TT|. Energetic tail distributions play a key role in coronal plasma 
dynamics |T2| and high resolution plasma observations near 1 AU confirm that even the 
distribution function of heavy solar wind ions are well fitted by a kappa-distribution ||13|| . 

The generation of velocity space distributions exhibiting pronounced energetic tails 
was frequently interpreted as a consequence of several different acceleration mechanisms. 
Those include besides DC parallel electric fields or field-aligned potential drops in re- 
connection regions also wave-particle interaction due to kinetic Alfven wave turbulence 
lis] and cyclotron interactions [|16|. Clear evidence of the importance of suprather- 



mal particle populations in space plasmas have motivated the development of a modified 



plasma dispersion function for kappa-distributions fLT\ , 

The use of the family of kappa-distributions to model the observed non-thermal fea- 
tures of electron and ion structures was frequently criticized since a profound derivation 
in view of fundamental physics was not available. A classical analysis addressed to this 
problem was performed by Hasegawa [|19] demonstrating that kappa-distributions turn 
out as consequence of the presence of suprathermal radiation fields in plasmas and Col- 
lier PU| considers the generation of kappa-like distributions using velocity space Levy 



flights. Furthermore, a justification for the formation of power-law distributions in space 
plasmas due to electron acceleration by whistler mode waves was proposed by Ma and 



Summers |21, 22, 03] and a kinetic theory was developed showing that kappa-like veloc- 



ity space distributions are a particular thermodynamic equilibrium state Here we 

demonstrate that the family of kappa-distributions results as consequence of the entropy 
generalization in nonextensive statistics, providing thus the missing link for the use of 
kappa-distributions from fundamental physics. 
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2 Theory 



A generalization of the Boltzmann-Gibbs-Shannon entropy formula for statistical equilib- 
rium was recognized to be required for systems subject to spatial or temporal long-range 
interactions making their behavior nonextensive. Any extensive formalism fails whenever 
a physical system includes long-range forces or long-range memory. In particular, this 
situation is usually found in astrophysical environments and plasma physics where, for 
example, the range of interactions is comparable to the size of the system considered. 
A generalized entropy is required to possess the usual properties of positivity, equiprob- 
ability, concavity and irreversibility but suitably extending the standard additivity to 
nonextensivity. 

In view of the difficulties arising in this conjunction within the Boltzmann-Gibbs 



standard statistical mechanics motivated Tsallis [^, ^ to introduce a thermo-statistical 
theory based on the generalized entropy of the form 



= kB—^ (1) 



1 - Ep I 

q-1 

where pi is the probability of the i*^ microstate, ks is Boltzmann's constant and q is 
a parameter quantifying the degree of nonextensivity and is commonly referred to as the 
entropic index. A crucial property of this entropy is the pseudoadditivity such that 

S,iA + B) = S,iA) + SgiB) + (1 - q)S,iA)S,iB) (2) 

for given subsystems A and B in the sense of factorizability of the microstate probabil- 
ities. For g 1 the standard Boltzmann-Gibbs-Shannon extensive entropy is recovered 
as 

Sg = -keJ^Pi^^Pi (3) 

Applying the transformation 

-^ = -n (4) 
q - 1 

to equation (|1]) and restricting to values — 1 < g < 1 yields the generalized entropy of 
the form 

= >^ks{Y.p]'"^ - 1) (5) 

for I < K < oo. With regard to Silva |^ one finds the corresponding one dimensional 
equilibrium velocity space distribution in kappa notation as 



f{v) = A, 



1 



(6) 
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where the normahzation constant reads 



A. 



N I V[n] 

vth v^r [k - 1/2] 



(7) 



Here denotes the particle density, Vth = \jlkBT jm is the thermal velocity where T 
and m are the temperature and the mass, respectively, of the species considered. The case 
K = oo corresponds to g = 1 wherefrom the Maxwell equilibrium distribution is recovered. 
Figure la illuminates schematically the non-thermal behavior of the distribution function 
(6) for some values of the spectral index kappa. 

Consistently, the one-dimensional equilibrium velocity space distribution in the q- 
nonextensive framework is written as 



1/(9-1) 



where 



Vth 



1-g- 



r[i/(i-g)] 



(9) 



r [1/(1 -g)- 1/2] 

for — 1 < g < 1. In analogy, the isotropic three dimensional velocity space distribution 
is found as 



/(V) = 5k 



1 + 



1 
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th 



(10) 



identical to the distribution function for a plasma in a suprathermal radiation field 
discussed by Hasegewa (1985) where the normalization constant reads 



5. 



N 



3/2^3^ /,3/2r[fi;_ 3/2] 
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TT 



for I < K < OO. Equation (|^) and (|T0D denote a reduced form of the standard kappa- 



distribution used for astrophysical applications and written conventionally as 



/(v) = A, 



1 



with the normalization constant 



A^ 



N 



7^3/2^3^ ^3/2 r [ft; - 1/2] 



(12) 



(13) 



For a variation of the spectral index kappa Figure 1 demonstrates the different behavior 
between the one dimensional case of the distribution function (12), generally used for 
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space applications and modeling, as compared to equation (6) resulting from the modified 
entropy approach. Since (6) and (12) differ only in the exponent no exact mapping is 
possible when substituting k, —>■ k, + 1. Nevertheless, highly similar structures are found 
for both with the appropriate values of kappa, hence favoring relation (6) for astrophysical 
velocity space distribution modeling due to the physical background. 



a) MODIFIED ENTROPY 




(b) TRADITIONAL KAPPA 
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Figure 1: Schematic plot of the family of kappa-distributions. The subplots (a) and 
(b) demonstrate the one-dimensional functional dependence of equation (6) and (12), 
respectively. For both, k = 2 corresponds to the outermost curve followed by values 
fi; = 3,4, 6 and 10. The innermost curve represents with k = oo an isotropic Maxwellian. 



An effective thermal speed 6 = Vthy {i^ — 5 /2) / k is commonly defind from the moments 
of the distribution function (12) where the spectral index kappa is conventionally limited 
to positive values k > 3/2 and we note that negative values of kappa have never been 
considered for space plasma applications. Contrary, the generalization to negative values 
of kappa identical to g > 1, discussed in the framework of nonextensive statistics as well, 
generates according to equation (^ or (|T^) a thermal cutoff at the maximum allowed 
velocity 



2kBT 



' K- 



m 



(14) 



hence providing an additional physical interpretation of the spectral index kappa of 
non-thermal plasma structures. To our knowledge there is presently no observational 
evidence for distributions corresponding to g > 1 in space plasmas available. On the 
other hand it should be mentioned that any observed thermal cutoff in velocity distribu- 
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tions would not have been recognized or interpreted as being a consequence of a possible 
nonextensivity of the system. 



Finally we note that Collier ||2^ has evaluated the Boltzmann entropy for kappa- 
distribution functions. Upon introducing the corresponding nonextensive velocity space 
distribution the resulting kappa-dependent entropy is found to obey qualitatively the 
same functional dependence. Since the suprathermal tails are more pronounced in the 
modified entropy CclSG clS compared to the standard kappa function, see Figure 1, also 
the entropy enhancement for low values of the spectral index, corresponding to a large 
fraction of suprathermal particle populations, is increased relatively to the traditional 
kappa-distribution case. Both approaches merge for large values of kappa towards the 
Maxwellian limit. 



3 Conclusions 

Nonextensive statistics was successfully applied to a number of astrophysical and cosmo- 
logical scenarios. Those include stellar polytropes |2^, the solar neutrino problem ||3U|| , 



peculiar velocity distributions of galaxies and generally systems with long range in- 



teractions and fractal like space-times. Cosmological implications were discussed ||3^ and 
recently an analysis of plasma oscillations in a collisionless thermal plasma was provided 
from q-statistics [^. On the other hand, kappa-distributions are highly favored in any 



kind of space plasma modeling |I| among others, where a reasonable physical background 
was not apparent. A comprehensive discussion of kappa distributions in view of exper- 
imentally favored non-thermal tail formations is provided by Leubner and Schupfer 



where also typical values of the index k, are quoted and referenced for different space 
plasma environments. 

In the present analysis the missing link to fundamental physics is provided within the 
framework of an entropy modification consistent with nonextensive statistics. The family 
of kappa distributions are obtained from the positive definite part | < k < oo, corre- 
sponding to — 1 < g < 1 of the general statistical formalism where in analogy the spectral 
index kappa is a measure of the degree of nonextensivity. Since the main theorems of 
the standard Maxwell-Boltzmann statistics admit profound generalizations within nonex- 
tensive statistics ||3^, ^ |3^, a justification for the use of kappa-distributions in 



astrophysical plasma modeling is provided from fundamental physics. 
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